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• . 1 INTRODUCTION 



ABSTRACT 

GRS 1915+105 has the largest mass function of any Galactic black hole system, al- 
though the error is relatively large. Here we present spectroscopic analysis of medium- 
resolution IR VLT archival data of GRS 1915+105 in the K-b&nd. We find an updated 
ephemeris, and report on attempts to improve the mass function by a refinement of 
the radial velocity estimate. We show that the spectra are significantly affected by 
the presence of phase-dependent CO bandhcad emission, possibly originating from 
the accretion disc: we discuss the impact this has on efforts to better constrain the 
black hole mass. We report on a possible way to measure the radial velocity utilising 
apparent _ff-band atomic absorption features and also discuss the general uncertainty 
of the system parameters of this well-studied object. 
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Since its discovery by ICastro-Tirado et alj (|l994l ) 
GRS 1915+105 remains one of the most intensively 
studied of all the Galactic X-ray sources, and is the 
prototypical Galactic 'microquasar'. This system is of 
particular importance, not only for our understanding of 
the formation of jets near black holes, but also because 
it is reg arded as a Galactic analogu e to Active Galactic 
Nucl e i iuMirabel fc Rodriguez! 1 1994 ICastro-Tirado et all 
1 19961 : iGreiner et al.1 l2001al lbTh Due to its location only 0.2 
arcsec from the Galactic plane, GRS 1915+105 suffers from 
strong interstellar absorption: the F-band extinction (A v ) 
is 19.6 ± 1.7 mag (|Chapuis fc Corbe]U2004n . For this reason 
observations in the if-band have been partic ularly critical, 
leadin g to the identification of the secondary l|Greiner et al.l 
l2001al). a radial ve locity estimate for the secondary 
( Greiner et al.l l2001bT) and the possib le determination of a 
photometric period (|Neil et al.l 120071 ). The radial velocity 
study yielded a mass function of 9.5 ± 3 .0 which , when 
combined with incl ination estimates (I Fender et al.l 1 19991 : 
IGreiner et afll2001bl ) indicated a primary mass of 14.0 ± 4.0 
Mq. This not only confirmed the presence of a black hole 
in the system but indicated it to be one the most massive 
black holes discovered to date in a Galactic X-ray binary. 
Other examples of massive primaries in clude Cygnus X-l 
with a primary mass of 14.8 ± 1.0 M (|Orosz et alj|201lh 
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12 + 2 Mq (IShahbaz et al.1 Il994h but has now been 



reduced to 9.0 + g;g Mq (|Khargharia et alj|2010h 



Foll o wing the estimate by IGreiner et al.l (|2001bl ). 
ICasaresI <l2007h noted that this mass falls outside the 
the oretical distribution rang e for our Galaxy as calculated 
by iFrver fc Kalogeral (12001 ). However, since the error on 



the black hole mass is so large (±4.0 Mq), it is clear that 
a more accurate estimate is required before we can make a 
meaningful c omparison to the theoret ical mass distribution. 
In addition, iMcClintock et all (120061 ) derive a lower limit 
for the dimensionless spin parameter of a, > 0.98 making 
GRS 1915+105 a near maximally spinning Kerr black hole. 
This conclusion was drawn from an analysis of RXTE 
and ASCA data for the thermal state of GRS 1915+105, 
and a model of the X-ray contin uum of a fully rel a tivisti c 
accretion disc. As discussed by IMcClintock et aT, I (|2006h . 
determining the mass of, and distance to, GRS 1915+105 
is of fundamental importance for verifying this method of 
spin determination. 



Here we present an analysis of archival Very Large Tele- 
scope (VLT) data an d a re-analysis of the data presented by 
IGreiner et al.|[2001al lbl In particular, we show in what fol- 
lows that, near orbital phase 0.25, the absorption lines are 
corrupted by emission, which will have a significant effect 
on the amplitude of the radial velocity variation. We discuss 
possible origins for this emission, and the effect it might have 
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Table 1. Complete list of archival data used including total exposure times for each night from 1999-2002. 



Date 




No. 


of Exp 


Total 


Phase 


Date 






No. of Exp 


Total 


Phase 


CUT) 








00 






(UT) 










(s) 




1999 


Jul 23 


10 x 300s 


3000 





46271 


2002 


Jul 30 


10 x 


300s 


3000 


0.23063 


2000 


Apr 23 


8 x 


240s 


1920 





74271 


2002 


Aug 


06 


20 x 


300s 


6000 


0.43595 


2000 


May 11 


8 x 


240s 


1920 





27410 


2002 


Aug 


08 


10 x 


300s 


3000 


0.49377 


2000 


May 23 


8 x 


240s 


1920 





63179 


2002 


Aug 


10 


10 x 


300s 


3000 


0.55101 


2000 


Jun 10 


8 x 


250s 


2000 





16171 


2002 


Aug 


12 


10 x 


300s 


3000 


0.63898 


2000 


Jun 13 


8 x 


250s 


2000 





25325 


2002 


Aug 


13 


10 x 


300s 


3000 


0.64005 


2000 


Jun 18 


8 x 


250s 


2000 





39898 


2002 


Aug 


14 


20 x 


300s 


6000 


0.65954 


2000 


Jul 03 


8 x 


250s 


2000 





83817 


2002 


Aug 


17 


20 x 


300s 


6000 


0.75728 


2000 


Jul 09 


8 x 


250s 


2000 





01589 


2002 


Aug 


19 


20 x 


300s 


6000 


0.81673 


2000 


Jul 12 


8 x 


250s 


2000 





10582 


2002 


Aug 


29 


20 x 


300s 


6000 


0.11109 


2000 


Jul 14 


8 x 


250s 


2000 





16342 


2002 


Aug 


31 


20 x 


300s 


6000 


0.17236 


2000 


Jul 17 


8 x 


250s 


2000 





25120 


2002 


Sep 02 


10 x 


300s 


3000 


0.22992 


2000 


Jul 23 


8 x 


250s 


2000 





46333 


2002 


Sep 05 


10 x 


300s 


3000 


0.28863 


2000 


Jul 27 


8 x 


250s 


2000 





54704 
















2000 


Aug 01 


8 x 


250s 


2000 





69426 
















2000 


Aug 21 


8 x 


250s 


2000 





28594 

















on the mass function and primary mass. We also discuss the 
feasibility of using //-band spectroscopy to measure the ra- 
dial velocity of the secondary from narrower atomic features 
using new Gemini observations. 

2 OBSERVATIONS 

2.1 ESO VLT archival data 

The European Southern Observatory (ESO) archival data 
used here were taken with the Infrared Spectrometer And 
Array Camera (ISAAC) on the 8-m VLT Antu telescope 
at Cerro Paranal (Chile). The short wavelength arm 
of ISAAC equipped with a 1024 x 1024 pixel Rockwell 
HgCdTe array with an image scale of 0.147 arcsec pixel -1 . 
Using the medium resolution grating (1.2 A pixel -1 in 
the if -band) yields a spectral resolution of ~3000 with a 
1 arcsec slit. Science data consisted of several 250-300 s 
individual exposures each night (see Table [T] for summary 
of the data used). We present a re-analysis of 16 nights 
of obser vations from 2000 M ay to August, which were the 
focus of[G rcmcr et al.ll2001al lbl and 13 more from 2002 July 
to September which are presented here for the first time. 

The bias subtraction, flat-fielding and sky subtraction 
were performed using the iraiQ CCDPROC package. The 
spectra were extracted using the iraf kpnoslit package. 
Initial wavelength calibration was carried out using arc 
spectra obtained before or after each science run. How- 
ever, cross-correlation of these spectra with the OH sky 
lines indicated a residual wavelength offset of ~3-5 A, 
presumably due to flexure in the telescope/spectrograph 
between the times when science exposures and calibration 
arc frames were taken. This wavelength shift was corrected 
for before the telluric features were removed. To correct for 
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Figure 1. Radial velocity curve of the previousl y unpublished 
2002 d ataset folded using the best-fit period of iGreiner et al.l 
ll2001bh P orb = 33.5 ± 1.5 days (marked with open circles). 
The fit overlaid has a reduced \ 2 = 1- The peak occurs at 
K d = 131. 5 ± 2.5 km s" 1 with a systemic velocity of 33.0 ±2.1 km 
s -1 . As can be clearly seen the p hasing of the dataset with the 
ephemeris of lGreiner et al ] l|2001ah leads to an offset in the peak 
from the expected 0.25 to ~0.5. Also plotted is the same data 
folded on our refined ephemeris (marked with crosses) from the 
combined dataset (sec Figure 3). Note that the refined ephemeris 
exhibits no significant offset from the expected peak of 0.25. 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 



telluric absorption the AO V star HD 179913 was observed 
either before or after the scie nce exposures each night 
and t he method as outlined by IVacca. Cushing fc Ravneri 
utilized. The data were exported to MOLIA0 where 
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Figure 2. (top left) L omb-Scargle peri odogram of the combined 2000 and 2002 dataset with a Tq = 2451666.5. (top right) Bretthorst's 
Bayesian periodogram l|Bretthorstl200lf l which has been marginalized to exclude noise resulting from a non-sinusoidal source, (bottom left) 
Normalized Log plot of the PDF to emphasise any lower peaks, (bottom right) Closer look at the two resulting major peaks which occur 
at 33.8 and 35.38 days respectively, with each having a FWHM of 0.1 days which we take as the error. Of these only the 33.8 day period 
successfully phases both datasets correctly. 



they were re-binned onto a common velocity scale. Three 
techniques were then used to measure the radial velocity 
curve: each spectrum of GRS 1915+105 was cross correlated 
against (i) a template spectrum of the Kill standard star 
HR 8117, (ii) the average of the (remaining) spectra of 
GRS 1915+105 and (hi) the average of the spectra occurring 
around phase 0.75 (see below). All three techniques resulted 
in similar velocities within the errors. 

As a hnal check, the spectra were also re-extracted using 
the VLT ISAAC pipeline, and the resulting radial velocities 
were again found to be consistent with those determined 
above once the shift due to instrument flexure had been 
removed. 

2.2 Gemini-North NIFS data 



In[G remer et all (|2001al ) the presence of atomic absorption 
features are identified in the r7-band. Motivated by the hope 
that a radial velocity measurement based on these inher- 
ently narrower features would lead to a more accurate and 
reliable mass function estimate, a total of three 'pathfinder' 
spectra were obtained using Gemini-North's Near-Infrared 
Integral Field Spectrometer on the nights of 2010 June 
29th, July 6th and July 13th. The observations were 
scheduled to sample the radial velocity variation at three 
orbital phases. NIFS provides 3D imaging spectroscopy 
with a 3.0" x 3.0" field of view, equipped with a 2048 x 2048 
Rockwell HAWAII-2RG HgCdTe array with an image scale 
of 0.103" x 0.043"per pixel across and along the slit respec- 
tively. Using the H-G5604 grating and the HK-60603 filter 



http: / / www. Warwick, ac.uk/go/trmarsh/ software 



yields a spectral resolution of ~5000. Science data consisted 
of 5 x 180s target exposures each night resulting in a S/N 
~24. The data were reduced using the iraf Gemini nifs 
package and telluric absorption corrected for by observing 
the A0V standard HD 182761 before and after science runs. 
The data were again exported to molly where they were 
re-binned onto a common velocity scale. The spectra were 
then cross correlated against the KI III template star HD 
83240, taken from the catalogue o f high resolution VLT 
IR spectra of iLebzelter et all ( 20121 ). A stringent mask of 
the spectra was used in each case to utilise only the lines 
identified in Figure [7] for the cross correlation. 



3 RESULTS 

3.1 ESO VLT archival data 

In Figure Q] we plot the previously unpublis hed data from 
2002 o nly, folded using the ephemeris of iGreiner et alj 
(|2001bl ). The maximum velocity shift is seen to occur 
at phase 0.5, whereas it should occur at phase 0.25 by 
definition, suggesting that a refinement in the ephemeris of 
IGreiner et aTT (|2001bl ) is needed. Both the 2000 and 2002 
datasets yield radial velocity amplit udes consistent with th e 
result of 140 ± 15 km s -1 found by IGreiner et al.l (|2001bl) . 
although the systemic velocity (7) differs somewhat to 
the previously published value of —3 ± 10 km s . This 



difference is most likely due to the skyline correction 
outlined in Section 2.1 as prior to this correction we 
find a systemic velo city that is consistent with that of 
IGreiner et all (|2001bl ) for the same dataset. 
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Figure 3. Radial velocity curve of the combined 2000 (marked 
with open squares) and 2002 (marked with crosses) datasets 
folded over the best-fit period of P or b = 33.8 ±0.1 days obtained 
from the Lomb-Scargle periodogram. The overlaid fit has a re- 
duced x 2 = !• The peak occurs at K<j = 131.9 ± 3.3 km s — 1 with 
a systemic velocity of 34.2 ±2.5 km s . The scaling of error bars 
to produce a reduced x 2 = 1 makes the 2002 error bars larger 
than those in Figure 1. This is due to the larger scatter of the 
points in the combined dataset. 



We next combined both the 2000 and 2002 datasets in 
order to derive an improved ephemeris by performing a pe- 
riod analysis via a Lomb-Scargle periodogram and then re- 
fining the results using Bretthorst's bayesian periodogram 
l|Bretthorstl 1200 ll ) to clean out any false signals, the re- 
sults of which are shown in Figure [2] ( calculated using 
T = 2451666.5 from lGreiner et al'l l|2001bl )V The peak that 
corresponds to a period of 33.8 ± 0.1 is the only one to suc- 
cessfully phase all the data consistently and so we adopt this 
as our revised orbital peri od for GRS 1915±105. This period 
is consistent with that of iGreiner et alj |2001bl ) (33.5 ± 1.5 
days) to within 1-er, but does not agree with the low am - 
plitude K-b&nd modulation discussed bv lNeil et al] (|2007h . 
Folding all the data on this revised period yields the radial 
velocity curve shown in Figure [3] A sinewave fit to this plot 
yields a secondary radial velocity semi-amplitude (Kd) of 
131.9 ± 3.3 km s _1 and a systemic velocity (7) of 34.2 ± 2.5 
km s _1 . To estimate the uncertainties on derived quanti- 
ties, we used a Monte Carlo analysis. This involved select- 
ing 100000 random samples from distributions of each input 
parameter i, Kd, Porb and Md (based on their means and 
standard deviations), and calculating the results based on 
the output distributions. Combining our initial estimates in 
P or b and Kd with the suggested donor mass (Md) of 0.8±0.5 



Mq (|Harlaftis fe Greiner |2004[) and th e binary inclination 



(i) of 66 ± 2 degrees jFender et aHll999h the mass function: 

frivn = ( M ° sin 0* = p orbK d 3 

{ ' ~ (M c + M D ) 2 ~ 2ttG 
is found to be 8.0 ± 0.6 Mq ( at the 68 per cent confidence 



interval), whic h is slightly lower th an the 9.5 ± 3.0 Mq 
determined by I Greiner et al] (|2001bl ). but consistent within 
their la error. The mass of the primary is found to be 
12 ± 1.4 Mq. 

However, in Figure [4] we plot the trailed spectra for the 
combined dataset phased on our refined ephemeris: emission 
can be seen in the first half of the orbital range, indicated 
by the lighter regions, blueward of all the CO bandhead 
absorption features (this is more apparent in the electronic 
colour version). The emission rises and falls over the first 
half of the orbit and is predominantly strongest near phase 
0.25. To demonstrate this emission more clearly, in Figures 
[5] & [6] we plot the summed spectra, phased on the revised 
ephemeris, for both the 2000 and 2002 datasets. In each case 
the data is further separated into two distinct orbital phase 
windows (0.15-0.35 and 0.65-0.85) to illustrate the presence 
and absence of the emission near 0.25 and 0.75 respectively. 
As can be clearly seen, there is significant CO bandhead 
emission present in the first half of the orbital phase range, 
not in the latter. This emission is present in the two epochs 
of observation (most noticeably near the 12 CO (2-0) and (3- 
1) transitions) indicating that it is not just a transient phe- 
nomenon. Whereas variability in the b andhead absorption 
had b een reported by previous authors (|Harlaftis fc Greinerl 
120041) . the fact that it occurs predominantly near orbital 
phase 0.25 has important implications for the radial veloc- 
ity curve, which we discuss in Section 4. 

3.2 Gemini-North NIFS data 

From our NIFS spectra (see Figure [7} we con firm the pres- 
ence of many of the atomic features noted bv lGreiner et all 
(|2001al ): namely the Al triplet ( A16723.5 A, A16755.2 A, 
A16767.9 A), Si I (A15964.4 A, A16685.3 A), Mg I (A15029.1 
A, A15044.3 A, A15770.1 A, A17113.3 A), Fe I (A15297.3 A, 
A15969.1 A, A16047.1 A) and OH (A16753.8 A). Equivalent 
widths of some of the major lines can be found in Table [2] 
We do not detect Br series (A or 77), CO or He I. We first 
used these spectra to independently estimate the effective 
temperature (T e g) of the secondary, from Equation 3 of 
iLe et al] (|201ll ). using the EW of the Mg I (A17113.3 A) 
line. From this we estimate an T e fi of 4540I5Q5 K which is 
reassuringl y consistent with the K 2 III classification already 
ascribed bv lGreiner et al] (|2001al ). 

The resulting velocities from our cross correlation anal- 
ysis (see Section 2.2 for details) are shown in Tabled These 
were then searched for the expected variation due to the 
motion of the secondary, based on the results of the it-band 
observations outlined in Section 3.1. However, even with the 
improvement in the error of the ephemeris (Section 3.1), over 
the intervening 110 cycles between the K & if-band obser- 
vations the uncertainty in phase information accumulates to 
~ ±0.33. Hence, in constraining Kd, the absolute phase was 
left as a free parameter. We used the 7 velocity found in Sec- 
tion 3.1. To our surprise, we do not find the radial velocity 
variation expected, based on the C O band-head measure - 
ments discussed in Section 3.1 and IGreiner et aT] (|2001bf ). 
Our measurements are consistent only with an upper limit 
of ~50 km s _1 on Kd, although this depends strongly on the 
assumed value of 7. 
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Wavelength (A) 

Figure 4. Trailed spectrum of all data, repeated twice for clarity. Lighter colours indicate emission and darker absorption. Notice the 
emission is visible only in the the first half of the orbital cycle, blueward of each bandhead, and is strongest near phase ~0.25 (see 
electronic version also). Intermittent emission at A23160 A is a residual telluric feature. 



Spectrum 


Mg (15770.9) 

(A) 


Si (1668.53) 

(A) 


Mg (17113.08) 

(A) 


Al (16723.5) 

(A) 


Velocity 
(km s _1 ) 


2010 Jun 29 
2010 Jul 06 
2010 Jul 13 
HD 83240 


2.02 ±0.16 
2.47 ±0.15 
2.55 ±0.15 
2.23 ±0.20 


0.59 ±0.10 
0.44 ±0.10 
0.49 ±0.10 
0.47 ±0.10 


1.07 ±0.10 
0.92 ±0.10 
1.05 ±0.11 
0.99 ±0.09 


1.25 ±0.22 
1.31 ±0.15 
1.07 ±0.13 
0.567 ±0.07 


9.1 ±20.5 

54.4 ± 21.3 

62.5 ±21.9 



Table 2. Equivalent widths of if-band spectra. Only the deepest lines are listed. Also listed are the corresponding radial velocity 
measurements, and their associated errors, from the cross correlation analysis in MOLLY. 



4 DISCUSSION 

H-Band Atomic Absorption 

It is clear that some caution should be used when interpret- 
ing the radial velocity variation (or lack thereof) in the H- 
band, as we have only three measurements at our disposal. If 
the lack of variation (compared to what is expected from the 
K-haxiA) is real, then we must look for alternati ve sources of 
the _ff -band emission. For ex ample, the work of lMuno et al.1 
|l999l ) and iLee et ail (|2002l ) suggests that GRS 1915+105 
may possess a circumbinary disc of material: if we assume 
such a disc is the source of the //-band absorption features, 
then we can use the FWHM of the latter to constrain the 
location of the emission. The observed FWHM of the ab- 
sorption features is ~60 km s , and assuming a Keplerian 
disc orbiting a binary of total mass 12.8 Mq, we estimate 
a radius of 530 R.0, well outside the binary radius of 100 
R,0. However, it remains suspicious that the i/-band spec- 
trum, despite the lack of radial velocity variation, so closely 
matches the expected spectral type of the secondary. It is 
clear that more systematic observations of these spectral fea- 
tures are warranted. 

CO Bandhead Emission 

The CO emission is a rare phenomenon in accreting 
binary systems: to our knowledge, it has only previously 
been observed in the well known cataclysmi c variable 
WZ Sagittae (iHowell, Harrison fc Szkodvl 12004 ). As with 
iHowell. Harrison fc Szkodvl (|2004 ). we also believe that H 
Pfund emission is unlikely to be a major contributor to 



the observed emission. Higher order lines (> P 33) blend 
to form a continuum at shorter wavelengths (see Figure [5]l 
and do not exhibit the discrete band structure blueward 
of each CO bandhead that is observed. Also, the intensity 
of Pfund emission lines should increase from shorter to 
longer wavelengths, whereas the emission that we observe 
is approximately uniform at each site. Hence we believe 
that CO bandhead emission is the best explanation for the 
observed emission features. 

With a dissociation energy of 11.1 eV, the site of this 
emission must reside in a relatively cool, outer part of the 
disc. Such a region cannot be directly illuminated by X-rays 
from the inner disc, and must instead be shielded from it. 
This could be provided, for example, by a warped accretion 
disc, which in turn has been invoked to explain the long 
term periodicitie s in the X-ray light curve of GRS 1915+105 
(|Rau et al.ll2003l ). 

However, it is the fact that the emission occurs prefer- 
entially near orbital phase 0.25 that most interests us here. 
Such emission may become more visible near this phase 
because the absorption lines of the secondary are at their 
maximum redshift at this phase, making emission from the 
accretion disc easier to observe. However, the same should 
be true at maximum blue shift (phase 0.75), which is not 
the case for our data. 

Alternatively, it may simply be that the shielded re- 
gion of the disc is best visible near binary phase 0.25 - 
although, in the context of a warped and precessing disc 
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Figure 5. Bottom: Average of 2000 spectra between 0.15 and 
0.35 in phase. Top: Average of 2000 spectra taken between 0.65 
and 0.85. Notice the emission present in the former compared 
to the latter. The positions of the Pfund lines are indicated to 
demonstrate that they cannot account for the structure of the 
observed emission. 
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Figure 6. Bottom: Average of 2002 spectra between 0.15 and 
0.35 in phase. Top: Average of 2002 spectra taken between 0.65 
and 0.85. Again, as with the 2000 spectra, the emission is present 
in the same orbital ranges. Again, the positions of the Pfund lines 
are indicated for reference. 

ijLee et al.ll2002r i. it is not clear why this should be the case 
for both epochs of observations. Another possible cause for 
the appearance of the blueshifted emission phenomenon is 
P-Cygni type emission originating fro m a wind driven out- 
flow of material, such as proposed by I Lee et all (|2002l ) for 
GRS 1915+105. This would imply a shell of ma terial some- 
wher e between 200 - 400 Rq from the primary (I Clark et al.l 
1 19991 . Equation 1). The dominant source of uncertainty for 
this estimate, however, is due to the uncertainty in the dis- 
tance to the binary. 

The uncertainty in the distance to GRS 1915+105 also 
plays a significant role in any attempt to tightly constrain 
the black hole mass. This is because the estimate for the 
primary mass stated in Section 3.1 relies heavily on the incli- 
nation of the system. While cited by many authors as 66 + 2° 
llMcClintock et al.|[2006l : Ivan Oers et~al1l2010| ; iRahoui et all 
l2010f ). this estimate is based in turn on an accurately 



known distance. However, the distance to GRS 1915+105 is 
still a matter of some deb ate; the conservative estim ate of 
9.0 + 3.0 kpc as proposed bv lChapuis fc Corbell (120041 ) would 
suggest a binary inclinatio n of 58 ± 11° (using Equation 

4 from iFender et alj Il999h. Combining this with a donor 
mass of 0.8 ± 0.5 M g (|Harlaftis fc Greiner1l2004h . yields a 
considerably more uncertain primary mass of 16.7+7.4 Mq. 

To estimate t he distance to GRS 1915+105, 
iGreiner et al] (|2001ri l combine d their estimate of 7 with 
the Galactic rotation curve of iFich. Blitz, fc Starkl (|l989h 
to yield D = 12.1 ± 0.8 kpc. Using the same technique, our 
estimate of 7 (see Section 3.1) yields D = 9.4 ± 0.2 kpc, 
where the low error is due to the fact that we assume a 
flat rotation curve and thereby only the uncertainty on 7 
contributes. This in turn would s u ggest an inclination of 
62 ± 2° (again from IFender et al.l (|l999h ), and a primary 
mass of 12.9 ± 2.4 M . 

Finally, it is also clear that the black hole mass esti- 
mate will also be affected by the presence of the bandhead 
emission. Specifically, near phase 0.25, the absorption line 
measurements will be skewed to higher velocities, artificially 
increasing the amplitude of the radial velocity curve. In ad- 
dition, the presence of the bandhead emission will also affect 
the measurement of 7 in the same manner. Hence, in this 
sense, the mass function we have derived must be considered 
an upper limit only. Further attempts to refine the mass 
function and 7 of GRS 1915+105 utilizing K-b&nd spec- 
troscopy will depend on the feasibility of decontaminating 
the absorption features from this emission contribution, es- 
pecially near phase 0.25. This will require higher resolution 
and S/N spectra than have been obtained thus far. 

5 CONCLUSIONS 

We have confirmed the pre sence of narrow atomi c features 
in the Zf-band as noted by IGreiner et al.l (|2001al ). but our 
analysis does not provide evidence for the expected radial 
velocity variation. This may suggest that the origin of the 
absorption is distinct from the secondary and outside of the 
binary but more observations are clearly required to con- 
firm this. Alternatively, the CO bandhead emission observed 
specifically at phase 0.25 brings into question the reliability 
of the radial velocity amplitude of the secondary based on 
the if-band observations. Higher resolution Zf-band obser- 
vations are required to disentangle the emission from the 
bandhead absorption so that a reliable radial velocity esti- 
mate can be made. 

ACKNOWLEDGEMENTS 

We gratefully acknowledge the use of the MOLLY software 
package developed by T. Marsh. We would also like to 
thank the anonymous referee for useful suggestions leading 
to an improved paper. 

REFERENCES 

Bretthorst G. L., Bayesian Inference And Maximum En- 
tropy Methods In Science And Engineering: 20th Interna- 
tional Workshop. AIP Conference Proceedings, Volume 
568, pp. 241-245, 2001 



The mass of the black hole in GRS 1915+105 7 

T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



id 

> 




J I I I I I I I I I I I I I I I I I I I I I I I L 



1.5x10 4 1.55x10 4 1.6x10 4 1.65x10 4 1.7x10 4 1.75x10 4 

Wavelength (A) 

Figure 7. .ff-Band spectrum confirming the presence of atomic features. Top: A spectrum from a single nights' observation with Gemini. 
All features indicated are found in each of th e three spectra. Equi valent widths of selected lines can be found in Table [2] Bottom: 
Template star HD 83240 from the catalogue of lLebzelter et al.l J2012h 

Casares, J., Black Holes from Stars to Galaxies - Across 

the Range of Masses. Edited by V. Karas and G. Matt. 

Proceedings of IAU Symposium, 238, pg3-12, 2007 
Castro-Tirado A. J., Brandt S., Lund N., Lapshov I., Sun- 

yaev R. A., Shlyapnikov A. A., Guziy S., Pavlenko E.P., 

1994, ApJS, 92, 469 
Castro-Tirado, A. J., Geballe, T. R., Lund, N. 1996, ApJ, 

461, L99 

Chapuis C, Corbel S., 2004, A&A, 414, 659 
Clark, J. S., Steele, I. A., Fender, R. P., & Coe, M. J. 1999, 

A&A, 348, 888 
Fender, R. P. , Garrington, S. T. , McKay, D. J. , Muxlow, 

T. W. B. , Pooley, G. G., et al., 1999, MNRAS, 304, 865 
Fich M., Blitz L., Stark A. A., 1989, ApJ, 342, 272 
Fryer, C. L., Kalogera, V. 2001, ApJ, 554, 548 
Greiner, J., Cuby, J., McCaughrean, M. 2001b, Nature, 

414, 522 

Greiner, J., Cuby, J., McCaughrean, M., Castro-Tirado, A., 

Mennickent, R. 2001a, A&A, 373, L37 
Harlaftis, E. T.; Greiner, J., 2004, A&A, 414, L13 
Howell, S. B., Harrison, T. E., Szkody, 2004, ApJ, 602, L49 
Khargharia, J., Froning, C. S., & Robinson, E. L., 2010, 

ApJ, 716, 1105 
Lebzelter T., et al., 2012, A&A, 539, A109 
Le, H. A. N., Kang, W., Pak, S., Im M., Lee, J.E., Ho, L. 

C, Pyo, T. S., Jaffe, D. T., 2011, arXiv. larXrv:1108. 14991 
Lee, J. C, Reynolds, C. S., Remillard, R., Schulz, N. S., 

Blackman, E. G., Fabian, A. C. 2002, Apj, 567, 1102 
McClintock, J. E.; Shafee, R.; Narayan, R.; Remillard, R. 



A.; Davis, S. W.; Li, L. X., 2006, ApJ, 652, 518 
Mirabel I. F., Rodriguez L. F., 1994, Nature, 371, 46 
Muno, M. P., Morgan, E. H., Remillard, R. A. 1999, ApJ, 

527, 321 

Neil, E., T., Bailyn, D., Cobb, B., E., 2007, ApJ, 647, 409 
Orosz, J. A., McClintock, J. E., Aufdenberg, J. P., et al., 

2011, ApJ, 742, 84 
Rahoui, F., Chaty, S., Rodriguez, J., et al., 2010, ApJ, 715, 

1191 

Rau, A.; Greiner, J., McCollough, M. L., 2003, ApJ, 590, 
L37 

Shahbaz, T., Ringwald, F. A., Bunn, J. C, et al., 1994, 

MNRAS, 271, L10 
Vacca, W. D., Cushing, M. C, Rayner, J. T. 2003, PASP, 

115, 389 

van Oers, P., Markoff, S., Rahoui, F., et al.; 2010, MNRAS, 
409, 763 



